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ABSTRACT 

X-ray reflection features observed from tfie innermost regions of accretion 
disks in Active Galactic Nuclei (AGNs) allow important tests of accretion theory. 
In recent years it has been possible to use the Fe Ka line and reflection continuum 
to parametrize the ionization state of the irradiated inner accretion disk. Here, 
we collect 10 measurements of ^, the disk ionization parameter, from 8 AGNs 
with strong evidence for reflection from the inner accretion disk and good black 
hole mass estimates. We find strong statistical evidence (98.56% confidence) for 
a nearly linear correlation between ^ and the AGN Eddington ratio. Moreover, 
such a correlation is predicted by a simple application of a-disk accretion theory, 
albeit with a stronger dependence on the Eddington ratio. The theory shows 
that there will be intrinsic scatter to any correlation as a result of different 
black hole spins and radii of reflection. There are several possibilities to soften 
the predicted dependence on the Eddington ratio to allow a closer agreement 
with the observed correlation, but the current data does not allow for an unique 
explanation. The correlation can be used to estimate that MCG-6-30-15 should 
have a highly ionized inner accretion disk, which would imply a black hole spin of 
~ 0.8. Additional measurements of ^ from a larger sample of AGNs are needed 
to confirm the existence of this correlation, and will allow investigation of the 
accretion disk/corona interaction in the inner regions of accretion disks. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: 
active — galaxies: nuclei 



1. Introduction 



It has been decades since it was realized that active galactic nuclei (AGNs) m ust be pow- 
ered by accretion onto supermassive black holes ( iLynden-Belll Il969l : lReeslll984| ). Yet, while 
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our theoretical understandi ng of the underlying phys i cs of accretion disks has evolved tremen- 



Balbus fc Hawlevlll998l : iBalbus 



dously over that time (e.g..lNoyikov fc Thornelll973l: IShakura fc Sunyaevlll973l : |Pringlelll981 



2003 



Hirose et al. 



20091 ). there are few direct measurements 



of the properties of AGN accretion disks that can be used to test the theories. The difficulty 
in detecting AGN accretion disks observationally arises naturally: they emit most of their 
radiation in the ultraviolet and optical, which can be easily absorbed and reprocessed by 
surrounding material and must also be separated from the starlight from the surrounding 
galaxy. When estimates of the disk continuum have been observationally determined (typi- 
cally from the outer part of the disk) they indicate broad agreement with the concept that 
the accretion dis k is an optically thick, thermal emitter with a radially dependent tempera- 



ture profile (e.g.. iKishimoto et al.ll2005l . |2008[ ). However, the few probes of the accretion disk 



emission from close to the black hole show signs of significant discrepancies with the the - 



oretical expectatioE0 (e.g., Koratkar fc Blaes 1999 : Davis et al. 2007 : Morgan et al. 2010 ). 



Accretion disk theory also predicts specific changes in the fiow p roperties (e.g., density, tern - 
perature) as a function of black hole mass and accretion rate (IShakura fc SunyaevI Il973l ). 
but there are currently very few observational tests of the s e predictions of the theory (e.g., 
Bonning et alll2007l : IPavis et al. l l2007l : iMorgan et aPboid IPavis k Laorlboilh . 



A direct probe of the inner accretion disks of AGNs is through the X-ray emission. 
All AGNs emit a fraction of their bolometric en ergy as a hard X-ray power-law (e.g., 
Mushotzky et al.l 1 19931 : IVasudevan fc Fabian! 120071 ) . The X-rays are thought to originate 
i n a hot and tenuous magnetically dominated corona that sits above th e surface of the disk 
jGaleev et al.lll979l : Iflaardt fc Maraschil E99ll . Il993l : Iflaardt et al]ll994[ ). Studies of this X- 
ray emission by generations of X-ray observatories discovered that in many AGNs there are 
spectral features superimposed on the power-law th at could be easily explained as repro- 
cessin g of the X-rays in nearby cold, dense material (jPounds et al.lll990l : iNandra fc Pounds 
19941 ). The strongest of these 'X-ray refiection' features are the Fe Ka emission line at 
6.4 keV, the associated Fe K absorption edge at 7.1 keV, and an over all hardening of the 
spectrum at high energies du e to electron scattering o f the X-ray photon s (iLightman &: White 
19881 : iGuilbert fc Reeslll988l : iGeorge k Fabian! 1 1 99 ll : iMatt et al.lll99lh . As the inner accre- 
tion disk likely subtends a large solid angle as seen from the X-ray source, it may dominate 
the X-ray refiection spectrum, although, depending on the source, contributions can arise 
from distant material such as the broad line region, the dusty absorber or the outer accre- 



^This should be contrasted with the situation for accretions disks around stehar-mass black holes, whose 
continuum is emi tted at X-ray energies and is, in general, well matched by the theoretical predictions (e.g., 



Davis et al.l 120061 ) . However, stellar mass black holes also present several p henomenological problems ( e.g., 



QPOs) that are yet to be fully explained by accretion theory (see, e.g. 
references therein). 
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tion disk ( iNandral 120061 : IShu et al.ll2010l ). High signal-to- noise X-ray spectral observations 



of some AGNs find a broadened, asymmetric Fe Ka line. The shape of this line is strong 



Tanaka et al. 


1995; 


Fabian et al. 


2002; 


Miller 


2007) 



spectra in these AGNs can probe physical properties of the inner accretion disk that would 
otherwise be extremely difficult to investigate by consideration of the disk continuum. The 
highly penetrating X-ray emission can be studied for nearly all AGNs and with very little 
contamination from the host galaxy. The only difficulty is to identify and model the com- 
ponent of the reflection spectrum that arises from the inner disk, as this will often be a very 



weak and low contrast featu re in the total spectrum (INandra et al.l 120071 ; iBallantynd 12010 
de la Calle Perez et allboioh . 



Once identified, the reflection sp ectrum from the inner a ccretion disk can constrain the 



i ron abundance of the a ccreting gas (IBa. 



( iNayakshin et al.l |2000| ; iBallantyne et al. 



lantvne et al. 



20021 ). the vertical density structure 



200 Ibl . I2OO4J ). and its ionization state, which is 



the focus of the current paper. The spectral features imprinted on the refle ction spectrum 
are a strong function of the ionization stat e of the illuminated 



_gas 



[e.g., 



Ross fc Fabian 



I993I : IRoss et al.l Il999t iRoss fc Fabian! l2005l : iGarcia fc Kallmaru 120101 ). In particular, the 
Fe Ka line transitions from a rest energy of 6.4 keV to 6.7 keV and then to 6.97 keV as 
the iron ions become progressively more ionized. In addition, the Fe K absorption edge 
moves to higher energies and Comptonization broadens out the natural width of the spectral 
features. The ionization of lower Z elements al so tends to flatten out the reflection spectrum 
( IRoss &: Fabianll2005l ; iGarcia fc Kallmanll2010l ). All these effects are potentially observable in 
high quality X-ray spectral data. Indeed, evidence for ionization of the disk surface has been 



ident i fied in the X-ray spectra of several AGNs (e.g. , iBallantvne et al. 



2OO4J : iLonginotti et al.l 120071 ; iBrenneman fc Revnoldsll2009l ; iNardini et al.l l20 111 ). 



2001a; Fabian et al. 



The photoionization of the surface of the inner accretion disk by the X-ray emitting 
corona can be parametrized by an ionization parameter, = ^-KFx/n-a, where Fx is the 
incident X-ray flux and tt-h is the hydrogen number density of the disk surface. A measure- 
ment of ^ can therefore provide information on the density structure of the disk and/or the 
illuminating conditions provided by the corona. Models of coronal generation predict that 
the a mount of accretion energy dissipated in the corona may be a function of t he accretion 
rate Jstella fc Rosnerlll984j ; iMerloni fc Fabianll2002l iBlackman fc Pessahl l2009[ l . and obser- 
vations show that the fraction of the bolometric lumi nosity released by an AGN in the X-ray 
band decreases with increasing Eddington ratio (e.g.. lWang et al.ll2004j ; IVasudevan &: Fabian 
20071). These two results imply there should be a dependence of ^ on the accretion rate (cf., 
Ballantyne &: Rossll2002l ). which, if measured, can t est how the density structure and corona 
power varies with accretion rate. Earlier work by llnoue et al.l (120071 ) indicated that ^ in- 
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creases with accretion rate, but those authors were unable to fit the X-ray data with ionized 
reflection models to measure ^ at various Eddington ratios. 

In this paper we measure the dependence of C, on AGN accretion rate by compiling the 
best available results from the literature. The sources were selected to have clear evidence 
of reflection from the inner accretion disk and robust black hole mass estimates. The next 
section describes the source selection in detail and presents the results of our experiment 
which are then discussed in the context of accretion and coronal models in Section [31 Our 
conclusions our summarized in Section HI 



2. Source Selection and Results 



The literature was searched to flnd the best available observationally determined values 
of ^. Sources selected for analysis must pass two criteria: flrst, the central black hole 
mass must be est imated either fro m reverberation mapping, or, if the radius- luminosity 
relationship (e.g., iKaspi et al.ll2005l ) is used, a careful decomposition of the H/3 hne must 



be performed to obtain an estimate of the virial motion of the broad-line region (BLR) gas. 
Second, the X-ray spectrum of the object must have strong evidence of reflection from the 
inner accretion disk (within 10 Vg, where Vg = GM/c^ is the gravitational radius of a black 
hole with mass M) in the form of a broad Fe Ka line. AGNs where relativistic effects are 
not needed to describe the line proflle are not included for analysis. These two requirements 
severely restricts the number of possible AGNs that can be included in the sample, but the 
objects that are included will provide the cleanest test of any relationship between and the 
Eddington ratio with a minimal amount of observational scatter. Finally, the AGNs must 
have a measurement of ^ from fltting a relativistically blurred ionized reflection spectrum to 
the Fe Ka line and coritinuu m. For consistency, this ^ must be estimated from the 'reflionx' 
model JRoss fc Fabian 2005 ). After searching through the literature, it was found that 10 



keV 



observations of 8 different AGNs satisfied the above criteria^ (see Table [J) . The 2- 
bolometric corrections, kx, are found from the measurements of IVasudevan &: FabianI (|2009| ) 
or, for three of the objects, are est imated by using the observed relationship between kx and 
photon index (jZhou fc Zhad 120101 ). 



Figure d] plots log^ versus the Eddington ratio for the objects listed in Tabled! and 



^Objects with relativistic Fe Ka lines that failed to make the list because of an uncertain black hole 
mass include MCG-6-30-15, IRAS 13224-3809, IH 0707-495 and Swift J2127.4-t-5654. Radio-loud AGNs 
with reverberation mapped black hole masses such as 3C 120, 3C 273 and 3C 390.3 were also not included 
to eliminate confusion due to emission from the relativistic jet. 
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shows a clear correlation between the two parameters (the simple linear correlation coefficient 
is r = 0.80). The Spearman rank correlation coefficient is 0.685 and gives a t- value of 
2.658, which corresponds to a 98.56% confidence level for 8 degrees of freedom. Thus, 
the observed correlation has strong statistical significance. However, the detection of the 
correlation depends on only three objects with (Lboi/-Z^Edd) > 0.2, and, as seen in Tabled! 
those AGNs have only estimates of the bolometric correction. Removing those three objects 
from the sample eliminates the correlation (r = 0.04); it is therefore crucial to obtain better 
measurements of bolometric corrections for high accretion rate AGNs in order to confirm 
the existence of the relationship between ^ and Lboi/-Z^Edd- 

To estima t e the functional relationship between these quantities, we follow the advice 



of llsobe et al.l (Il990l ) who recommend the ordinary least-squares bisector method. This 
calculation finds the following, nearly linear, relationship between and the Eddington 
ratio: 

loge = (1.008 ± 0.162) log(Lboi/i:Edd) + (3.14 ± 0.164). (1) 

This fit is shown as the dashed line in Fig. [H and is formally consistent with the slope (0.805) 
and intercept (2.95) calculated from a simple least-squares fit. If the observed errors in ^ 
are included than this line gives a reduced » 1; this can be reduced to ~ 1 if all the 
errors in log^ are increased to ±0.4. However, as discussed below, substantial scatter in this 
relationship is likely to arise naturally, so a reduced ~ 1 would not be expected. 



3. Discussion 

3.1. Physical Interpretation 

Figure [1] clearly shows that AGNs that are accreting at larger fractions of their Edding- 
ton rate exhibit more ionized inner accretion disks. As mention in Sect. 1, such a relationship 
might be expected depending on how accretion energy is dissipated in the corona at different 
accretion rates. In this Section, we use simple a-disk accretion theory to explore the possible 
physical origins of this observed relationship. 

As the ^ measurements were all chosen to arise from the innermost regions of the accre- 
tion disk, radiation pressure will likely dominate the support of the accretion flow. Likewise, 
relativistic effects will also be non- negligible. Finally, it will also be necessary to include 
the effects of a non- zero (and potentially variable) fraction of accretion energy dissipated in 
the corona, /. The ionization parameter can be approximated as ^ ~ rripLx/H'^p, where 
rrip is the mass of a proton, Lx is the total X-ray luminosity incident on the accretion disk 
which has a gas density p, and H is the distance from the X-ray source in the corona to the 
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reflecting region of tlie disk. As tliere is no complete understanding of the structure of the 
corona close to the black hole, for simplicity we will assume a geometrically thick corona 
so that H/R = 1, whe r e R i s the radial distance along the disk. Theoretical arguments 



by iBlackman fc Pessahl (120091 ) show that the survival of a magnetically dominated corona 
will require large s cale magnetic fields, so a thick corona above the inner disk may not be 
unrealistic (see also lMerloni fc Fabianll200ll ). With this assumption. Appendix lAl shows that 
for a radiation pressure supported disk. 



^ - (4.33 X 10») (^) - (ii 



^bol 



^Edd 



R 



-7/2 



— R^'^Rj}R\f{l 



f)^ erg cm s 



(2) 



wher e r] is the radiative efficiency of the disk, a is the viscosity parameter (jShakura fc Sunyaev 
19731 ). and {Rr, Rz, Rt) encompasses the general relativistic effe cts and are dimensionles s 



functions of a*, the dimensionless black hole spin, and (R/rg) (e.g.. lNovikov fc Thorndl 19731 ). 
This estimation of ^ is nicely independent of the central black h ole mass. Initially, we assume 
constant values of a = 0.1 and t] = 0.089 (jPavis fc Laoiil201ll ). so a model curve of ^ versus 



(Lboi/-^^Edd) can be described by only three parameters: a*, R/rg and /. 

Figure [2] plots several models calculated from equation [2] and compares them against the 
data and the least-squares bisector fit from Fig. [H There are three sets of curves for three 
different (a*, R/rg) pairs. As described in detail below, the solid lines and the long dashed 
line differ in the treatment of /. Fig. [2] is remarkable in that the simple a-disk model of 
Eq. [2] predicts reasonable values of ^ for realistic values of the parameters; there is no reason 
for this to be expected a priori. This basic level of agreement provides confidence that Eq. [2] 
can be used to explore the physical reasons behind the observed correlation between ^ and 

{Lhol/ -^^Edd)- 

There can be a wide range of expected values of ^ for a given value of LboiZ-f'Edd de- 
pending on the value of the black hole spin and the radius which dominates the reflection 
signal. As these values will likely vary from object to object, and the latter may vary from 
observation to observation, a fair amount of natural scatter will occur in the observed rela- 
tionship. In fact, this scatter may be a significant contribution to the observed slope of the 
relationship; however, the scatter is not so significant as to destroy the correlation. This fact 
argues that AGNs within this range of Eddington ratios are not uniformly distributed over 
and i?, but are confined to a relatively small range in these parameters. Clearly, more 
ionized refiection fits in a wider range of AGNs is needed to better understand the intrinsic 
scatter in this relationship. 



Vasudevan &: Fabian! (120071 ) make use of their measured AGN X-ray bolometric correc- 
tions to estimate that the coronal fraction / ^ 0.45 for low Eddington rate AGNs. The 
long dashed line in Fig. |2] shows the predicted ^-(Lboi/-^^Edd) relationship from eq. [2] for this 
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value of / when a* = 0.998 and R = 8 rg. This hne has a slope of 3, as predicted from 



eg. El much steeper than the measured slope of ~ 1 (eq. [T]). However, IVasudevan fc Fabian 
( 120071 ) also found that / is lower for high Eddington ratio sources, dropping to ~ 0.11. The 
anti-correlation of / with Eddington ratio is a natural consequence of the well establishe d 
diminishing of relative X-ray power in high luminosity quasars (e.g., ISteffen et al.ll2006l ). 
Mod els of accretion disk corona predict different dependencies of / on the Eddington ratio, 
with IWang et al.l (12004! ) providing a comparison with observa tional estimates. The s teepest 
dependence of / with (i^boi/-^Edd) is provided by the model of IStella fc Rosnen ( 1l984l ) which 
finds / oc (LboiZ-^^Edd)' 
/ = 0.45 at (Lboi/^Edd 

still significantly steeper than the observed slope. 



^''^'^ . Inserting this dependence into eq. 121 and normalizing so that 
= 0.01, yields the solid lines in figure |2l which have slopes of ~ 2.3, 



Although the intrinsic scatter is likely to contribute to the mismatch between the ob- 
served and predicted slope, it is interesting to consider what other physical mechanisms 
may be operating that can bring the model closer to the observed correlation. There are 
numerous possibilities, but, as outlined below, not all are likely physical reasonable. 



3.1.1. Coronal Fraction 

If / oc (Lboi/-Z^Edd)~^ then the slope of the model is approximately equal to the observed 
one. However, this dependence is so strong that, if normalized to 0.45 at (Lboi/-^^Edd) = 0.01, 
/ drops to 1% at an Eddington ratio of only 0.07. Such a strong dependence of / would be 
very difficult to reconcile with the observed X-ray luminosities of luminous quasars. 



3.1.2. Gas Density 



Perhaps the simplest change to the model is to alter the gas density of the disk. The 
value of p used in calculating ^ is an approximate, vertically-averaged value of the gas density, 
and may not be entirely accurate for the outer layers of the disk that is being illuminated 
by the X-rays. Furthermore, the X-ra y heating may cause the disk scale height to inc rease 
and thus decreasing the density (e.g., iNayakshin et al.ll2000t iBallantyne et al.ll2001bl ). As 
there is more accretion power emitted as X-rays at low Eddington ratios, then its plausible 
that this density correction could be dependent on (Lboi/-Z^Edd)- The dotted line in figure 121 
shows the a* = 0, = 8 model after reducing the density with a factor that varies as 
(Lboi/-Z^Edd)~^'^- The reduction in density ranges by a factor of ~ lO'' at {L\^oi/ L^dd) = 0.01 
to 25 for an Eddington ratio of unity. 
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There are problems with this solution, however. First, it can most easily be applied 
to low spin black holes. Higher spin black holes, such as the a^, = 0.998 models plotted 
in figure |2] would require a denser disk at high Eddington ratios to bring the model into 
agreement with the data. In addition, refiection calculations of disks which do adjust into 
hydr ostatic balance and have low density surfaces predict that ionized Fe Ka lines are com- 
mon (IBallantyne fc Rossll2002l ) unless the disk is very weakly illuminated. It is not clear that 
this is consistent with the observed variation in coronal fractions, where more accretion power 
is dissipated in the corona at low Eddington ratios and the ionization parameter is observed 
to be small. Both these problems may be mitigated if processes within the disk (e.g., photon 



bubb l e or Parker instab ilities) transported denser material to the surface (IBallantyne et al. 



2OO4J : iBlaes et al.l 120071 ). More simulations investigating the density structure of accretion 



disk photospheres are needed to quantify the possible corrections to the a disk densities. 



3.1.3. Radiative Efficiency 



According to equation [21 ^ oc r^^^, so if the radiative efficiency of accretion disks was 
proportional (Lboi/-Z^Edd)''^^° then, after including the effects of the variable /, the slope of 
the model would be more in line with the observed correlation. Of course, 77 is a measure 
of the total binding energy available to be radiated away in the accret ion fiow and thus 
depends on the black hole spin and the innermost radius of the disk (e.g., lKroliklll999l ). The 



maximal values of t] are limited to be 0.31 (0.038) for maximally spinning bl ack holes that 
are co -rotating (counter-rotating) with respect to their accretion disks [see iDavis fc Laor 
( I2OIII ) and references therein]. Additional physics can alter these expectations, however; for 
example, magnetohydrodynmical effects may provide an extra torque to the gas within the 
traditional innermost stable or bit and increase the efficiency above the theoretical upper- 
limit (e.g.. lAgol &: Krolikll2000l ). Alternatively, accretion energy may be advected inwards or 



used to dri ve outfl o ws wh ich can reduce the efficiency. Observationally, arguments based on 
the work of ISoltanI (jl982[ ) have shown that the average efficiency of accreting black holes over 



their accretion history is > 0.1 (e.g.. Yu fc Tremainel 2002; Elvis et al. 2002; Marconi et al. 



2004; Barger et al. 2005 ). Recently, Davis fc Laoi ( 201ll ) estimated 77 for 80 PG quasars and 
found an average value of ~ 0.1 (with a wide scatter) and no strong correlation with the 
Eddington ratio. Indeed, implementing the rj cx (i^boi/-^^Edd)^^^° dependence into equation [2] 
finds that it can only hold over a factor of ~ 20 in Eddington ratio before reaching the 
theoretical limits. 



Interestingly, iDavis fc Laorl (j201l[ ) find a positive correlation between 77 and black hole 
mass, ?7 oc M^/^. Thus, an alternative way to reduce the slope of the ^-(Ivboi/-Z^Edd) relation- 
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ship is with (M/Mc^) oc ( L^o] / L^mY^^ ■ Again, such a strong relationship is not observed 
( lKollmeierll2006l : ISteinhardt fc Elvisll2010l : iDavis fc Laorll201ll ). nor is it even hinted at with 
the data in Table [T] where only a slight negative correlation is possible (r = —0.44). 



3.1.4- Disk Viscosity 



The models shown in fig. [2] assume a constant a 
relationship could be reduced t o ^ 1 if a oc (L hoi /-^Edd 



0.1, but the slope of the predicted 
Observational estimates of 
a typically find a > 0.1 (e.g.. King et al.l 120071 ). but num erical simulations of accretion 



fiows typically find values an order of magnitude lower (e.g.. IDavis et al.ll2010l ). The strong 
relationship between a and the Eddington ratio required to account for the observed ^- 
(Lboi/-^^Edd) relationship results in either a large fraction of low Eddington AGNs with a ~ 1, 
or that most rapidly accreting AGN have a < 0.01. As the a parameter is at best a 
parametrization of our ignorance of the accretion process, it is difficult to assess the validity 
of such a relationship between a and accretion rate without guidance from simulations that 
explore disk viscosity ClS db function of M. 



3.1.5. Black Hole Spin 

Figure [2] shows that the black hole spin has an important impact on the predicted values 
of ^, and that the slope of the predicted curves could be reduced if a* was a decreasing 
function of (Lboi/-^^Edd); an idea that has some theoretical justification. In particular, some 
models of black hole growth predict that only the largest black holes that have stopped 
rapidly accreting have developed high spins, while the lowest mass black holes which ma y 



still undergo rapid accretion will generally have low values of a* (e.g.. iFanidakis et al.ll201ll ). 
That is, the high Eddington ratio AGNs would be lower mass black holes with small spins, 
and the low Eddington ratio AGNs would be high mass black holes with large spins. This 
scenario i s also a possible exp lanation for the correlation between t] and black hole mass 
found by iDavis fc Laorl (120111 ). Although the objects in our sample do not present strong 
evidence for a correlation between black hole mass and Eddington ratio, we assumed that 
oc {L\,o\l L-EAdY searched for an x that could reduce the predicted slope to match the 
observations. The minimum slope found was 1.4 for x = —0.11 and a refiecting radius very 
close to the innermost stable circular orbit. Thus, a dependence of spin with Eddington ratio 
cannot, on its own, reproduce the observed correlation between ^ and (Lboi/-^^Edd)- 
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3.1.6. Other Options 



There are other physical effects neglected in equation [2] that could inffuence the observed 
^. F or example, if the X-ray source is close to the black hole, gravitational light bending 



(e.g., 



Miniutti &: Fabian! l2004f ) will increase the irradiation of the disk and thus increase ^. 
If this was more common, or more effective, at smaller (Lboi/-^^Edd) (as might be expected 
given the larger values of /), then this effect could reduce the slope of the model to be closer 
to the observed value. However, if the disk is more highly illuminated than this would also 
increase the strength of the reflection spectrum and increase the equivalent width (EW) of 
the broad Fe K g lines. Iron lines with large (i.e., several hund red eV) EWs are rare in 
Seyfert galaxies ( Nandra et al. 2007 : de la Calle Perez et al. 201o[) which argues that strong 
light bending is not a common phenomenon (see also iBallantynd l2010l ). Alternatively, it 
is possible some fraction of the accretion energy dissipated in the corona is used to drive 
an outflow (e.g., iMerloni fc Fabian! |2002[ ). This would have two important effects for the 



predicted value of first, this would reduce the energy available to illuminate the accretion 
disk, and, second, the X-rays emitting region may be moving rapidly away from or towar d the 
disk and thus significantly decrease or increase the irradiating flux ( iBeloborodovl 119991 ). As 
/ is likely to depend on Eddington ratio, then it is certainly plausible that such effects may 
also be stronger over a certain range of (Lboi/-^Edd). More generally, it is possible (perhaps 
likely) that all of the effects considered in this section are in play at some level in real AGNs. 
Measurements of ^ from the inner regions of accretion disk are needed from several more 
AGNs with reverberation mapped black hole masses to further explore the nature of this 
correlation. 



3.2. Application to MCG-6-30-15 

One potentially important application of this correlation is that it may help distinguish 
between two competing interpretations of an observed X-ray spectrum. For example, MCG- 



analyzed many times since 


the line was 


irst detected by ASCA f( 




Tanaka et al. 


1995; 


Ballantvne & Fabian 


2001 




Wilms et al. 


2001; 


Fabian et al. 


2002; 


Ballantvne et al. 


2003; 


Revnolds et al. 


2004; 


Brenneman & Revnoldsl 


2006; 


Miniutti et al. 


200/ 


When modeled 



as one powerful relativistic line, reflectio n fits to the spectra of MCG-6-30-15 yield log^ 5, 2 
and require a large iron abundance (e.g.. lBrenneman fc Reynoldsll2006l : iMiniutti et al.ll2007l ). 
However, equally good fits to the spectra can be found if the red wing of the line arises from 
highly ionized Fe close t o the black hole and the neu t ral 6.4 keV core arises f r om much larger 
distances on the disk ( Ballantyne fc Fabian! l200ll : iBallantyne et al.l l2003l : [Reynolds et al. 
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2004 ) . None of the AGNs in Table [T] show evidence for a strong second reflector from the 
accretion disk that could bias the measurement of the ionization parameter and the resulting 
correlation. Thus, this degeneracy is unlikely to be common, but, given the importance of 
MCG-6-30-15 to the field, it is interesting to use the correlation observed in Fig. [T]to provide 
additional information on the likely ionization state of the MCG-6-30-15 accretion disk. 

Unfortunately, the black hole mass of this sour ce is not well constrained with a current 
best estimate of 4 x 10^ Mq jMcHardv et allboosh . The flux of MCG-6-30-15 whe n it was 



observed in a deep minimum state (2.3 x 10 erg cm ^ s ^; IWilms et al.l 120011 ) is used 
to estimate a lower limit to log^. As there is no direct measurement oi kx, we utilize the 



photon index of MCG-6-30-15 in the deep minimum state (F = 1.8; iRevnolds et a. 



2004) 



and again make use of the correlation between kx and F described by IZhou &: Zhad (120101 ) 
which results in kx ~ 24. Combining these measurements results in Lboi/-^Edd = 0.18 which 
when inserted into Eq. [1], and taking the extremal values of the slope and intercept, yields 
log^ > 2.1. This lower limit suggests that the accretion disk of MCG-6-30-15 is most likely 
to be significantly ionized. 

Interestingly, when the broad Fe Ka line of MCG -6-30-15 is fit by a high l y ionized reflec 



tor, t he ionization parameter is usually log^ > 3.5 (IBallantyne et al.ll2003l : [Reynolds et al. 
20041) ■ and such a highly ionized line is naturally broader due to Comptonization (e.g. 



Ross et al.lll999l ). Thus, an ionized accretion disk in MCG-6-30-15 would reduce the value of 

g^. req uired to explain the width of the Fe Kg l ine (e .g. jReynolds et al.ll2004j : iBrenneman fc Reynolds 
20061 ) . However, the fits of IBallantyne et al.l (120031 ) found that the ionized reflector was con- 
strained to a radius of i? ~ 5 r^, and therefore a non-zero black hole spin will be necessary. 
We can be more quantitative by combining the measured ^ and R from spectral fits with the 
predicted ^-(i^boi/-^Edd) correlation. An ionization parameter of log^ > 3.5 would be larger 
than expected for the correlation seen in Fig. (H but scatter due to the various values of a^: 
and R will bring objects ab ove that line. Accord ing to Equation [21 if log^ > 3.5 at i? ~ 5 
for Eddington ratios ~ 1 (IMiniutti et al.l 120071 ). then a* ~ 0.8. Thus, a rapidly spinning 
black hole is still necessary for MCG-6-30-15, but a very large value is not required if the 
inner disk is ionized. 



4. Conclusions 

In this paper we have presented evidence from the literature that the ionization state of 
the inner accretion disk in AGNs is correlated with the Eddington ratio of the accretion flow. 
Very conservative criteria were employed to select the data from the literature; namely, there 
must be strong evidence for reflection from within 10 Vg, the reflection spectrum must have 
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been fit with the latest ionized disk models, and there must be a high quality estimate of the 
central black hole mass. Despite these steps, the correlation relies on estimates of the X-ray 
bolometric corrections for the high Eddington rate sources. Therefore, precise and accurate 
measurements of bolometric corrections (especially of rapidly accreting AGNs) are needed 
to both confirm the existence of this correlation and to better define its intrinsic scatter. 

The existence of a relationship between ^ and i^boi/-^Edd is predicted by simple a disk 
accretion theory, although with a slope that is steeper than what is observed. It is not 
possible to determine an unique explanation that can bring the two slopes into agreement, 
although a Eddington ratio dependent decrease in the disk density (possibly related to the 
changing coronal power) presents the simplest possibility. If the correlation can be better 
defined by future observations than it may lead to insight into the physics of the energy fiow 
between the disk and X-ray emitting corona. 

Finally, we showed that the ^-(LboiZ-^^Edd) correlation will also be useful in distinguishing 
between different interpretations of AGN X-ray spectra. As an example, given the current 
estimate of the black hole mass of MCG-6-30-15, the correlation predicts that the accretion 
disk of MCG-6-30-15 should be significantly ionized with log^ > 2.1. This result provides 
evidence that the strong line observed from this object may arise from two distinct reflecting 
regions with different ionization states, rather than one low-.^ area with a very high Fe 
abundance. The high values of C, measured in this interpretation require a* ~ 0.8, reducing 
the value of the spin necessary to explain the width of the red wing. 

We thank the referee. Dr. Chris Reynolds, for a helpful report that improved the paper. 
This work was supported in part by NSF award AST 1008067 to DRB. 



A. The Ionization Parameter in a Radiation Pressure Dominated Disk 



To estimate the ionization parameter of an accretion disk illuminated at a radius R by 
a X-ray luminosity Lx, we re- write ^ = AirFx/nn as ^ ~ nipLx/ H'^p, assume H/R = 1 (see 
Sect. 13. ip , and use the density of a radiation pressure dominated disk provided by (IKrolik 
1999h : 



p = (2.23 X 10 



-6\ 



0.1 



a 

oT 



1 / M 



Me 



-1 



L 



bol 



^Edd 



-2 /RX 3/2 

{^] RlRTR],^l-f)-'gcm-\ 



VgJ 



(Al) 

where the corr ection due to a non-zero coro nal dissipation fraction / has been included from 
the results of ISvennson fc Zdziarskil ( 119941 ). The relativistic corrections {Rr, Rz, Rt) are 
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simple, but lengthy, analytic functions of a^, the dimensionle ss black hole spin, and the disk 
radius (R/rg) and can be found elsewhere (e.g.. iKrolikI 119991 ) . 



Finally, it is assumed that that the X-ray luminosity is powered entirely by the fraction 
of the accretion energy that is dissipated in the corona; i.e., Lx = /-Z^boi- That is, none of the 
coronal energy is used to launch an outflow or jet, which is a reasonable as sumption given the 
relati vely high Eddington ratios of the objects making up the sample (e.g. 
2OO2I ). Combining this assumption with eq. lA II yields equation [2J 



Merloni fc Fabian 



B. The Ionization Parameter in a Gas Pressure Dominated Disk with 

Electron Scattering Opacity 

Although not likely applicable to the inner disk reflection observed in the AGNs con- 
sidered here, for completeness we also derive an estimate for ^ for the case of an illuminated 
ga s pressure do minated disk (with electron scattering opacity). In this case, the gas density 
is jKroliklll999l ): 

^ - <-) iur'° 

(Bl) 

where again the correction due to dissipation in the corona was taken from the work of 



Svennson fc Zdziarskil (jl994[ ). Then, with the same approximation as before for the produc- 



tion of the ionizing X-rays (i.e., Lx = /Lboi), the ionization parameter is: 

^ - ii^r i^r i^r it^T e. e. ... 

(B2) 

The high density of a gas pressure dominated disk predicts a very low ^ for supermassive 
black holes, and with a much slower dependence on the Eddington ratio than the radiation 
pressure dominated disk. 
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Fig. 1. — A plot of log^ versus LboiZ-f'Edd for the AGNs listed in Tabled A significant corre- 
lation is observed with a Spearman rank correlation coefficient of 0.685 (98.56% confidence 
level). The dashed line plots ordinary least-squares bisector fit to the data (Eq. [1]). 
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Fig. 2. — As in Fig. [H but overlaid with model curves derived from Eq. El To guide the eye, 
the gray dashed line show the least squares bisector fit from Eq. [T] and Fig. [T] The solid lines 
plot the predicted ^-(i^boi/-^Edd) relationship from Eq. [2]for different combinations of a*, the 
black hole spin, and R, the radius of p eak reflection. These particular models assume / oc 



(Lboi/IvEdd)"°-^^ (IStella fc R,osnedll9841 ). The long dashed line plots t he a* = 0.998, R = 8 



model if a constant / 
shows the a* = 0, i? = 



0.45 is assumed ( IVasudevan fc FabianI 120071 ). The dotted black line 



8 Vg model when the disk density has been reduced by a factor that 



is proportional to (Lboi/-^^Edd)~'^'^- This factor ranges from ~ 10^ for -Lboi/-^Edd = 0.01 to 25 
at i^boi/-^Edd = 1- As described in the text, this is only one out of several possibilities to 
bring the slope of the model into agreement with the observed relationship. 



Table 1. Details of AGN Sample 









s 


^ Ref. 


^ J 


Lx Ref. 








Mrh Ref 




Mrk 766 


NLSl 


0.0129 


061+329 


BR09 


7.6 X 10*2 


N07 


57 


ZZIO 


1.7 X 10^ 


WIO 


2 


NGC 3783 


Sy 1.5 


0.0097 


on+24 
'J'J-20 


BR09 


1.1 X 10*3 


N07 


15 


VF09 


2.8 X 10^ 


WIO 


0.04 


NGC 4051 


NLSl 


0.0023 




BR09 


2.0 X 10*1 


N07 


67 


VF09 


1.7 X 10^ 


DIO 


0.06 


Ark 120 


Sy 1 


0.0327 




BR09 


9.4 X 10''3 


N07 


25 


VF09 


1.4 X 10^ 


WIO 


0.13 


Ark 120 






280lg 


Nil 


7.5 X 1043 


Pll 










0.10 


Fairall 9 


Sy 1.2 


0.047 




Pll 


1.3 X 10*4 


Pll 


10.5 


VF09 


2.6 X 10^ 


P04 


0.04 


Mrk 79^ 


Sy 1.2 


0.022 


177+^ 


Gil 


1.7 X 10*3 


Gil 


10.5 


VF09 


5.3 X 10^ 


P04 


0.03 


Mrk 79*^ 






60^1 


Gil 


9.5 X 10*2 


Gil 










0.01 


I Zw 1 


NLSl 


0.0611 


306011°^ 


G07 


4.3 X 10*3 


G07 


76 


ZZIO 


2.0 X 10^ 


MIO 


1.3 


Mrk 478'=''^ 


NLSl 


0.079 


746^1™ 


Z08 


2.3 X 10^3 


G06 


86 


ZZIO 


3.2 X 10^ 


MIO 


0.48 



Note. — Details compiled from the literature of AGNs with high quality black hole masse estimates that also hav e X-ray spectra 
that e xhibit clear evidence for a relativistic Fe Ka line and have been fit with the 'reflionx' ionized disk model of iRoss fc Fabian 
Lx denotes the intrinsic rest-frame 2-10 keV luminosity of the AGN, Mbh is the mass of the black hole, and Lboi is the 
estimat ed bolometric luminosity . The Eddingt on luminosityj LEdd; is define d as -^Edd = 1 - 3 x 10 3^(MBH /;^f^) erg s~^. Ref erences: 

IWoo et all torn. VF09 



BR09: 



Brenneman fc Revnoldsl toO^. m7: iNandra et al.l hoojl. ZZIO: Izhou fc Zhaol J2010l). WIO 



Vasudevan fc FabianI J2009l). D loT TPennev et all |2010l). Pll: IPatrick et all (I2OI1I), N il: iNardini et al.l ( 2 01ll) P04: IPeterson et all 
tooi ). Gll: iGallo et al.l feoill) . G07: iGallo et al.l (l2007h . MIO: iMarziani et al.l toidt) . Z08: IZoehbi et al.l (l2008l) . G067[Gallo, (.2006 ). 



^data taken from the long Suzaku observation 



'^data taken from the XMM6 observation 



■^data taken from longest XMM-Newton observation 

iGallol I2OO6I ) listed the 2.5-10 keV flux, so Lx wfll be slightly underestimated 



